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Analytical and numerical arguments are presented in case of a pair of two-state systems in a singlet 
state that the threshold visibility for testing Bell's theorem on the entire range of measurement 
settings is 33.3 %. It is also shown that no lower treshold exists. 
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Bell's theorem states that local realistic theories can- 
not give account of all the quantum mechanical correla- 
tions. The best known proofs are based on the Clauser- 
Horne |Q or the Clauser-Horne-Shimony-Holt j^j inequal- 
ity. These inequalities are consequences of local realism, 
thus experimental demonstrations of their violation ||, 
j| imply the necessity of a revision of at least one fun- 
damental physical concept. 

This far reaching consequence motivates the continued 
search for new methods for formulating Bell's theorem 
and testing local realism. A recent direction is to exploit 
the full angular dependence of the quantum mechanical 
correlations, which is equivalent with the case of infinitely 
many apparatus settings at each side of the experiment 
|§- 0- A standard Bell type experiment is assumed, 
i.e., a source emits two spin-half particles in a maximally 
entangled (singlet) state, so that the particles do not in- 
teract after the emission. Subsequently, the projection of 
the spin at a direction a is measured on the first particle 
and the projection of the spin at a direction b is measured 
on the second particle (a, b are unit vectors). In a real 
experiment one has losses and thus the joint probability 
of having the results m and m' (— ±1) is 



Pqm(iti, m'; a, b) = — (1 — mm'Vab) 



(1) 



Here < V < 1 stands for the visibility. It is an im- 
portant question which visibility suffices to falsify the 
premises of local realism in case of infinitely many ap- 
paratus settings. In a recent paper Zukowski showed 
that the threshold visibility is not larger than 3/4. This 
improved the previous estimates 8/7r 2 || (derived for the 
case of coplanar settings) and ir/4 ||. In the present pa- 
per we argue that the threshold visibility is actually 1/3, 
and no lower value is possible. Throughout we adhere to 
the notations of Ref. . 

Local realism implies that the above joint probability 
is of the form 

P HV (m,m';a,b) = / dXp(X)P A (m\a, X)P B {m , \b, X) (2) 

J A 

where A stands for the 'hidden variable' characterizing 
the common past of the particles and A for the space of 
its allowed values. The question is whether the actually 



observed joint probability (|l]) admits the representation 
(^) . We shall present an analytical and a numerical argu- 
ment. The analytical treatment gives a precise value for 
the threshold visibility, but involves some additional as- 
sumptions concerning the nature of the hidden variable. 
In case of the numerical treatment we do not rely upon 
these additional assumptions, but the resulting value of 
the threshold visibility is determined only up to a finite 
accuracy. The analytical result coincides with the nu- 
merical one within this accuracy. 

In case of the analytical argument we make the follow- 
ing three assumptions: 



1. As we have identical particles, 

PA(m\n,X) = P B (-m\n, X) . 



(3) 



2. As PA(m\n, X) is a scalar, rotational invariance sug- 
gests that A is actually a unit vector, A and 



PA(jn\ri, A) = f(mriX) 



(4) 



3. Rotational invariance also requires that there is no 
distinguished direction, thus 



p(X)dX = —dtt , 
Air 



(5) 



where d£l stands for an infinitesimal spatial angle 
around the direction A. 

Let us represent the function f(x) (cf. Eq.(|]), note 
that — 1 < x < 1) as a series of Legendre polynomials: 



/(x)=^c J P J (x) (6) 
j=o 

Here Cj-s are real coefficients. 

Inserting Eqs.(|)-(|) into Eq.(|) and using the identity 



-> 4-7T 

Pj (aA)P fc (bX)dil - 5 jk ——Pj (3b) (7) 



we have 



1 



P H v{m, m'\ a,b) = ■ T Pj(-mm'ab) (8) 

3=0 ^ + 1 

On the other hand, Eq.(Q) can be written as 

PQAiim, m'\ a, b) = —P (—mm'db) + —P\(—mm!ab) (9) 

Comparing this with Eq. (Q) we get that cq and c\ are the 
only nonzero coefficients and 



a = ±- 



'3V 



(10) 
(11) 



If a and A are perpendicular, Pi{—mm!a\) is zero, thus 
the positivity of Pernio, A) requires (cf. Eqs.(||), (||)) 
that in Eq.(|l0|) the positive sign must be chosen. If 

A = ±a, then Pi(—mmaX) — ±1, thus the positivity 
of P A (m\a, A) requires that c ± ci > 0, i.e., by Eqs.(^0|), 

(O) 



^3 



(12) 



As |Pi(x)| < 1, Eq.(12) ensures the positivity of 
Pvl(m|a, A), Pb(?7i'|&, A) for any directions a, b, A. This 
readily implies that the threshold visibility is 1/3, as be- 
low this value the observed joint probability ([!]) admits a 
local realistic representation (||), while above this value 
the positivity of the conditional probabilities Pernio, A), 
Pb(to'|&, A) cannot be achieved. 

Although the assumptions (||)-(|5|) look reasonable, one 
might suspect that the above result hinges upon these 
assumptions and perhaps without them the joint proba- 
bility ([!]) admits a local realistic representation even at 
higher visibilities. We present a numerical argument that 
it is not the case. 

The numerical procedure begins with choosing N di- 
rections for a and N directions for b and then equating 
Eqs.(Q) and (||). In the latter the integral is estimated 
by a sum, i.e., the equation 

1 - M 

-(I - mm'Vajbk) = PnPA (m\aj , ti)Pb (m'\bk , n) (13) 

71=1 

is to be solved. This is equivalent with 

M 



where 



Vajb k = ^2 PnAj, n B k 



A j<n = 1-2 P A (l\Sj,n) 
B k , n = 2P B (l\b k ,n)-l 



(14) 



(15) 



\Aj,n\ < 1 

\B k J < 1 



(16) 



and 



^ PnAj, n = 

n=l 

M 

^ PnB k . n = . 



(17) 



At the solution of Eq.(|14|) one can utilize the singular 
value decomposition 



(18) 



i=l 



where Ujj and V kl i are orthogonal matrices. Note that 
there are only three nonzero singular values p n , owing to 
the three dimensionality of the vectors Sj , b k ■ A solution 
of Eqs.([l4|)-(17) can be then obtained in the following 
way: 

1. One chooses three M dimensional vectors qi and 
another three M dimensional vectors £j and makes 
them orthogonal to each other, 



and to the vector (y/pi, ^/~p~2, ■■■^/~p^) T , i.e. 

n 

^ \J~pl~i(U)n = . 



(19) 



(20) 



2. Calculate 



A 'j,n = U j,iy/Pi(%)n/VPn 
i=l 
3 



(21) 



3. A solution of Eq.(14) is given by 

A jtn = WA' j>n 
B Kn = Wb' 



4. Eq.(gJ) implies that 

l/V7 = max{|4,„|,|P/J} 



(22) 



(23) 



which satisfy 
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Applying this scheme, a Monte-Carlo simulation has been 
used to find the threshold visibility. The steps of the 
simulation are the following: 

1. Choose the unit vectors Sj, bk at random 

2. Choose the vectors ql, U at random and normalize 
them according to Eqs.(p^), (ptj|). 

3. Choose p n at random and normalize by 

5>™ = 1 ( 24 ) 



4. Calculate V from Eqs.(|T]), (H|. 

5. Change a component of <fi, ti or p n at random and 
repeat the calculation of V . If the resulting value 
is larger than the previous one, keep the changes, 
otherwise discard them. 

6. After having found the maximal V for fixed Sj, bk, 
change these vectors at random and find (by re- 
peating the previous steps) the corresponding V 
again. The minimum of these V values is selected, 
because at the corresponding setting for dj, bk a 
local realistic representation of the observed joint 
probabilities at higher visibility cannot be given. 

Obviously, one is interested in the limit N — > oo. For a 
finite N one has less restrictions, thus the resulting esti- 
mate for the threshold visibility is higher. By increasing 
TV the estimates decrease monotonously. As for the num- 
ber M, numerically it turned out that increasing its value 
makes the convergence slower but does not influence the 
results for V. Hence one may set it to the minimum 
M = 4. The values for N have been changed from 3 to 
10 4 . The convergence proved to be rather slow, e.g., the 
estimate for the threshold visibility at N = 1000 is still 
0.37 ± 0.001. The final numerical result (extrapolation 
for N — > oo) is 0.33 ± 0.03. This is consistent with the 
previous analytical result. 

It is instructive to compare these results with those one 
may obtain from Bell's inequality and the Clauser-Hornc- 
Shimony-Holt inequality. Since these inequalities do not 
exploit the full angular dependence, one gets a higher 
estimate for the threshold visibility. Note, however, that 
at the derivation of Bell's inequality 

P HV (l,l;a,b) + P HV {l,l;b,c) > P„ v (l,l;a,c) (25) 

strict anticorrelation (i.e., PhvO-i 1; b, b) = 0) implied by 
the singlet state is also assumed. Inserting Eq.(|l|) into 
Eq.(p5[) we have 



V 
~2 



(3- (a + c-6) 2 ) 



< 1 . 



(26) 



2/3 « 0.667, as below this value Eq.(|26|) is always satis- 
fied, while for a larger value of V it may be violated for 
suitably chosen directions a, b, c. This value is smaller 
than the previously known smallest value 3/4 = 0.75 [Q, 
but here the additional information about the proper- 
ties of the singlet state has also been utilized. This is 
equivalent with the assumption that in the special case 
of a = b (cf. EqM) the visibility is 100%. Using only 
the requirements of local realism together with the full 
angular dependence we obtained above the much lower 
threshold visibility l/3~ 0.333. 

Let us consider now the Clauser-Horne-Shimony-Holt 
inequality 

P HV {1, 1; a, b) - P HV {1, 1; a, b') + P HV (1, 1; a', b) 
+P H v(lA;ct,b') - Phv(1;Z') - P HV (l;b) <0 (27) 

Here P H v{^,a') = J2 m > p hv{1, m'; a', b) and 
P HV&b) = Y, m p Hv{m,l;a',b). Inserting Eq.(Q) into 
Eq.([27|), the resulting inequality can be cast to the form 



^ ((a + b' - b) 2 + {a' - b' - b) 2 - 6) < 2 



(28) 



For fixed V, b and b' the l.h.s of inequality ( |28| ) is the 
largest if a is parallel to b' — b and a! is antiparallel to 
b' + b. Then the lengths of a and b' — b add up to the 
length of a+b'~b (similarly in case of a'). Thus, denoting 
the angle between b and b' by tp, we get from Eq.(^8|) the 
inequality 



2V2Vi 



- + - < 2 

2 4/~ 



(29) 



This readily implies that the l.h.s. is maximal (for a fixed 
V) if if = 5 and the threshold visibility is not larger than 
l/\/2 « 0.707. Note that even this value is smaller than 
the previously obtained smallest value 3/4 = 0.75 0. 
The CHSH inequality relies only upon the assumption 
of local realism, but does not exploit the full angular 
dependence. As we have seen, taking into account the 
latter, too, a much lower value for the treshold visibility 
can be deduced. 
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The l.h.s. will be the largest when a + c — b = 0. This 
implies that the threshold visibility V is not larger than 
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